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We analyze the transport through a narrow ballistic superconductor-normal-metal-super- 
conductor Josephson contact with non-ideal transmission at the superconductor-normal-metal in- 
terfaces, e.g., due to insulating layers, effective mass steps, or band misfits (SIN interfaces). The 
electronic spectrum in the normal wire is determined through the combination of Andreev- and 
normal reflection at the SIN interfaces. Strong normal scattering at the SIN interfaces introduces 
electron- and hole-like resonances in the normal region which show up in the quasi-particle spec- 
trum. These resonances have strong implications for the critical supercurrent I c which we find to be 
determined by the lowest quasi-particle level: tuning the potential /j, x o to the points where electron- 
and hole-like resonances cross, we find sharp peaks in I c , resulting in a transitor effect. We compare 
the performance of this Resonant Josephson- Transistor (RJT) with that of a Superconducting Single 
Electron Transistor (SSET). 

PACS numbers: 74.80.Fp, 74.50. +r, 74.60.Jg 
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I. INTRODUCTION 



The ability to control the supercurrent flow through 
narrow superconductor-normal-metal-superconductor 
(SNS) contacts is not only of scientific interest-but also 
provides many opportunities for applicationsa. With 
recent progress in nanofabrication technology it has be- 
come possible to study devices in which electrons prop- 
agate ballistically and where the transport proceeds via 
few conduction channels!! Using gated SNS junctionsH, 
the transparency of the normal region can be manipu- 
lated: in a transparent wire, Andreev scattering at the 
NS boundaries produces phase-sensitive quasi-particle 
levels which carry large supercurrents; conversely, if the 
transmission is not ideal the admixture of normal scat- 
tering reduces the supercurrent transport. Accordingly, 
SNS junctions with a tunable transmission through the 
normal part defiue_a natural setup for a superconducting 
transistor deviceatl. 

Recent interest on transport through narrow chan- 
nels and quantum point contacts concentrates on di- 
verse phenomena suck m conductance quantization in 
normal constrictionsa>EliL2l, supercpfr^Bt quantization in 
superconducting SNS junctionsB&EJilij, or the transistor 
effect in super-links using .either gated structuresaErt-3 
or injection tcchniqueslljllZl. Experimentally, such 
junctions are fabricated uajng—superconductor — sexni- 
conductor-heterostrucfprje^iJoEfl, break-junctionsLTLl, 
metal-nanplithographyE3llZl, or with the use of carbon 
nanotubesEEI. Theoretically, transport through no rmal 
constrictions has been studied by Glazman et al£3t3 
within a quasi-classical description assuming adiabatic 
joints between the channel and the leads. The corre- 
sponding extension to superconducting leads by means 



of a scattering matrix approach is due-, to Beenakker0, 
see also Ref. EOL while Furusaki et alS3 proceeded with 
the numerical analysis of junctions with non-adiabatic 
geometries and non-ideal interfaces. The evolution of 
the quasi-particle spectrum and the supercurrent quan- 
tization in a gated narrow SNS junction, as well as its 
transformation into a SIS tunnel junction- has been re- 
cently described by Chtchelkatchev et al&3 Here, we ex- 
tend this analysis to the study of SINIS junctions, where 
T stands for a non-ideal interface between the supercon- 
ducting banks S and the normal channel N. 

In an ideal SNS junction the quasi-partiele spectrum 
is determined by the Andreev-scatteringE2l at the SN 
boundaries producing phase-sensitive levels transporting 
large supercurrentstS. The position and relative arrange- 
ment of these states strongly depends on the chemical 
potential in the wire as well as on the phase difference 
between the superconducting banks. The inclusion of 
weak normal scattering at the SIN interface will only 
softly modify this quasi-particle spectrum through the 
mixing of (left- and right-) current carrying Andreev 
states. On the contrary, strong interface scattering in- 
troduces electron- and hole- like resonances within the 
normal region, defining a new starting point. The posi- 
tion of these resonances again depends on the wire's ef- 
fective chemical potential. Tuning a pair of electron- and 
hole-like resonances to degeneracy these will be mixed by 
the Andreev-scattering and new phase sensitive levels are 
formed carrying supercurrent. This mechanism then pro- 
vides a natural setup for the implementatioB_of a super- 
conducting Josephson field effect transistoiU'Q where the 
supercurrent is switched on and off by tuning the scatter- 
ing resonances into degeneracy through the manipulation 
of a potential, e.g., through external gates. Model stud- 
ies of such systems for individual resonances have been 
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carried out recently by Wendin and coworker; 




FIG. 1. Narrow channel SINIS contact with non- ideal in- 
terfaces, e.g., insulating layers, effective mass steps, or band 
misfits. We consider adiabatic constrictions to avoid the mix- 
ing of transverse channels, (a) Geometrical setup showing 
the gates narrowing the wire, (b) potential landscape with a 
flat barrier bounded by 5-scatterers modelling insulating lay- 
ers. All sources of non-ideal transmission account for normal 
reflections (NR), whereas the discontinuities of the gap pa- 
rameter A are responsible for Andreev reflections (AR). 



with Ho = —h 2 d x [l/2m(x)]d x + U(x) — [i x (x) and where 
u v and v u denote the electron- and hole-like components 
of the wave function The potential U(x) induces 

normal scattering at the NS interface and is due to an 
insulating layer or a band offset, for example. Simi- 
larly, the mass m{x) may change at the NS interface, 
again generating normal scattering. The transverse en- 
ergy e±(x) of the— channel is enclosed in the effective 
chemical potentially fi x (x) = e F — £j_(a;). We will also 
make use of the kinetic energies E = e F ± e of electron 
(+) and hole (— ) states as measured with respect to the 
band bottom in the superconductors. 

The spectrum splits into continuous and discrete con- 
tributions and we will concentrate on the latter part with 
e v < A in the following, as it provides the main contribu- 
tion to the critical supercurrent in the most interesting 
transport regimes, see below. We solve the Bogoliubov- 
de Gennes equation for the SINIS junction with the help 
of the usual transfer-matrix techniques. In order to do 
so we first have to determine the resonance structure of 
the related NININ problem that arises if we replace the 
superconductors by normal metallic leads. 



A. NININ-junctions 



Below we proceed in two steps: After a brief definition 
of the problem (Sec. II), we first determine the scatter- 
ing states of the corresponding problem where the super- 
conducting banks are replaced by normal- metallic leads, 
the NININ junction. Second, we reinstall the supercon- 
ducting banks and determine the mixing of the (electron- 
and hole- like) normal scattering states through Andreev- 
scattering at the normal-superconductor interface. In 
Sec. Ill we analyze the quasi-particle spectrum for weak 
and strong normal scattering at the NIS interfaces. Sec- 
tion IV is devoted to the calculation of the supercurrent; 
we discuss the functionality of the Resonant Josephson 
Transitor (RJT) and compare this device with the Su- 
perconducting Single Electron Transistor (SSET). 



We expand the scattering states in a basis of in- and 



outgoing states \p™ ° ut on both si< 
phases that vanish at the interfaces 



of the wire, with 
±L/2, 



* L (x) = a£ e *(*+£/2) +a o« e -<fc(x+L/2) ) 



tf H (a;) = ane- lk(x - L/2) 



(2) 
(3) 



withk(E) = V2mE/h the wave vector of an incident par- 
ticle. The energy dependent scattering matrix S connects 
the expansion coefficients a^ ° ut (here, t and r denote the 
moduli of the matrix elements), 



*exp(«x*) 
rexp(ix r ) 



-rexp(i[2x* - x r ]) 
iexp(ix*) 





< 







(4) 



II. BALLISTIC CONTACTS: SCATTERING 
MATRIX APPROACH 

We consider a narrow metallic lead with few transverse 
channels connecting two superconducting banks with a 
rectangular pair potential A(x < —L/2) = Aexp(i<p L ), 
A(|x| < L/2) = and A(x > L/2) = Aexp(i^ R ), see 
Fig. 1(a). Joining the channel adiabatically to the su- 
perconducting-banks, the transverse channels in the wire 
are separableEJa and for each of them the quasi-particle 
spectrum e v is determined through the ID Bogoliubov- 
de Gennes equation (we choose states with e v > 0) 



(1) 
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with this definition of the 5-matrix and the basis 
states, the scattering phases account for the propagation 
through the wire. E.g., for a vanishing transverse energy 
in the normal region and in the absence of any interface 
barrier a finite phase x = kL is picked up. 

We describe the effective chemical potential in the 
normal region through a smooth function characterized 
by its minimum fi x (0) = [i x q and a positive curvature 
mf2 2 = d x /i x (see Fig. 1(b)). The parameter fi x o (i.e., 
the diameter of the metallic wire) is assumed to be tun- 
able by means of external electrostatic gates. In a long 
wire the potential is flat, hfl <ti A, and produces a sharp 
switching between transmission and reflection within the 
energy interval Tifi; the corresponding transverse energy 
s±(x) defines a smooth and flat potential barrier in the 
interval (—L/2, L/2). In the following we refer to such 
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barriers as smooth and tiexp(ix ) and exp(ix ri ) de- 
note the associated 'inner' scattering amplitudes describ- 
ing the motion of the quasi-particles between — L/2 + 
and L/2— (we use an analoguous definition of the 'inner' 
scattering amplitudes as in (H)). The particular smooth 
geometry of such inn(p,barriers justifies the application 
of the Kemble formulalij for the transmission probability, 
tf = l/{l + exp[—27r(^ x (0)±£)/?ifi]}, whereas the quasi- 
classic method can be used to determine the scattering 
phases x ti an d X Ti ■ 




-10 
E- e ± (0) [meV] 

FIG. 2. Modulus and phase of the transmission amplitude 
through a parabolic barrier with <5-scatterers at the edges (pa- 
rameters: Z = 2, £j_(0) = 5 eV, L — 2500 nm). Note that the 
energy E is swept and £x(0) is kept fixed. A similar resonance 
structure is found for other sources of non-ideal transmission 
at ±L/2. 

In order to obtain the global scattering amplitudes 
iexp(ix') and rexp(ix r ) of the junction we have to in- 
clude effects of normal scattering at the boundaries. As 
a typical example we consider symmetric smooth barri- 
ers of height e F — [i x0 bounded by Dirac scatterers of 
strength Vq at the interfaces (such a setup describes an 
NTNTN junction with insulating layers at the interfaces; 
see the Appendix for a brief discussion of other typical 
cases). 

First, we concentrate on the global transmission am- 
plitude iexp(ix') which is easily expressed in terms of 
the scattering amplitudes of the inner barrier and the 
Dirac <5-scatterers, see the Appendix for details of the 
calculation. The main effect of the scattering at the 
channel boundaries is to introduce resonances at ener- 
gies E rcs which we call perfect (imperfect) if t(E rcs ) = 1 
(< I). For energies E > £j^(0) + HQ/2 the electrons 
easily propagate through the channel, while the reflec- 
tion at the boundaries produces almost equidistant reso- 
nances with a spacing A w ir/dEX ti ~ ^(p-xo^^) 1 ^ 2 , with 
e L = h 2 TT 2 /2mL 2 (see Fig. 2; the last relation describes 
the case of a flat inner barrier). On the contrary, for 



small energies below e±(0) — h~£l/2, tunneling suppresses 
the propagation through the channel and only imperfect 
resonances separated by ~ 2A survive. In the interme- 
diate region, pairs of perfect resonances collapse and be- 
come imperfect. Moreover, the scattering phase x grows 
with the energy and picks up a phase of 7r and 27r at per- 
fect and imperfect resonances, respectively. We refer to 
resonances at energies E ICS = £f + £res above e F as elec- 
tronic and to those below, at E les = e F — e ros , as hole-like 
resonances. Below, we will be mostly concentrating on 
the propagating regime with E > £j_(0) + Ml/2. In this 
regime the global transmission amplitude takes the form 
(U± « 1) 



iexp(ix*) = 



I - Z 2 + i2Z + Z 2 e i2 * 



(5) 



with Z = mVo/h fc F the dimensionless parameter giving 



the strength of the scattering potential, 
width 



Resonances of 



r = 



9eX 1 



2A T 
~2-T 



with 



A = 



d E X u 



(6) 



and a transmission probability T = 1/(1 + Z 2 ) appear as 
the denominator in (|^) touches the complex unit circle 
(note that T — ► 2h~v/L, with v the particle velocity, in 
the absence of any scattering potential) . 

For symmetric barriers the reflection amplitude 
rexp(ix r ) is determined (up to ±1) by the unitarity of 
the scattering matrix, i.e., r = (1 — t 2 ) 1 / 2 and x r 
\ ' + Tt(n + 1 /2). The integer n jumps by unity at perfect 
resonances. 



B. SINIS-junctions 

After evaluating the normal scattering amplitudes we 
can reinstall the superconductors and match the scat- 
tering states in the normal region with the evanescent 
modes in the superconducting banks. We make use of 
the Andreev ap , pr |Qxi mat ionE3 and obtain the quantiza- 
tion conditiontiltj'CJ 



cos(x+ — X- ~ a ) = r+r- cos (3 + i + t_ cos ip, 



(7) 



where the +(— ) signs refer to the-.kinetic energies e F ± e 
of the electron(hole)-like statesE3. The Andreev scat- 
tering at the NS boundaries introduces the phase a = 
2arccos(e/A) decreasing from tt at e = to at the 
gap £ = A, as well as the phase difference ip = <^ L — ip R 
between the two superconducting banks. For symmet- 
ric barriers the phase [3 = (x\. — X+) ~ (x~ — X-) i s a 
multiple of tt and produces a smooth function r + r_ cos (3 
changing sign at perfect resonances [see Fig. 5(b)]. 

The case of an ideal SNS junction where the ef- 
fective chemical potential joins smoothly to the band 
bottom in the supepeonductors has been analyzed by 
Chtchelkatchev et alx3 In this special situation the global 
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and inner scattering amplitudes coincide and resonances 
are absent. Sweeping fi x o, they find a one parametric 
(fJ, x o) family of discrete spectra describing the transition 
of an insulating SIS tunnel junction (at fj, x o <C —A) into 
a ballistic SNS structure (at fi x o » A), see Fig. 3(a). 
More precisely, electronic levels with an exponentially 
small dependence on the phase ip are converted into 
phase sensitive Andreev-levels as /i x q is increased, i.e., 
as the band bottom in the wire is lowered. It turns 
out that the critical supercurrent is carried by the low- 
est state, I c = ma,x tp [(2e/h)d ip e ], and is realized for a 
phase tp = 7r — 0. Increasing the channel width, the 
critical supercurrent increases in steps of e/(ro + ?i/A) 
as new transverse channels open. The travel time tq of 
the quasi-particlcs is easily calculated within the quasi- 
classical scheme and approaches the asymptotic value 
To ~ L/v F . x in the open channel. This analysis explains 
the dependence of the non-universal critical supercurrent 
steps on the junction parameters. 

In the following we go beyond the analysis of 
Chtchelkatchev et al. and determine the spectrum and 
the supercurrent transport in a narrow ballistic SINIS 
Josephson junction including resonant barriers in the 
normal region, see Fig. 1(b). 



III. QUASI-PARTICLE SPECTRUM 

Below we shall see that the bound state spectrum of 
junctions with smooth barriers is only slightly modified 
by switching on weak resonances, see Fig. 3(b), with 
the degeneracy of the Andreev levels at ip = and 
7r lifted. For a fixed pair of levels this splitting 

fc = e' + ) — is roughly a periodic function of /u x o 
with a period A reflecting the resonance spacing. As the 
strength of the interface scattering increases the discrete 
spectrum is gradually distorted (see Fig. 3(c)); the phase 
sensitivity at large chemical potential fj, x o (region I) is 
reduced and the bound state energies align with the res- 
onance energies of the normal NININ setup. The effect 
of interface scattering on the continuum part of the spec- 
trum has been analyzed in Ref. E3l 



A. Weak Resonances 

We first concentrate on weak resonances characterized 
by a small scattering parameter Z <C 1. In particular, 
we aim at estimating the level splitting in the Andreev 
spectrum shown in Fig. 3(b). Using Eq. (|^) we evalu- 
ate the terms on the right side of the quantization con- 
dition (pi) to lowest order in Z and find = 1 — 
2Z 2 [cos^x+ + cos 2 X-] and r + r_ = 4Z 2 1 cos x+ 1 1 cos x'_! | • 
For symmetric barriers and as a consequence of the 
unitarity of the scattering matrix, the rescaled phase 
(3(e) /it is an integer valued function starting from 
at e = and jumping by unity whenever e lines up 



with a (perfect) resonance at e rcs or e los of the junc- 
tion. We then may write cos/3 = (—1)™, where n(e) de- 
notes the number of perfect resonances within the inter- 
val (0,e). Hence, cos/3 generates a smooth contribution 
r + r_ cos/3 = 4Z 2 cos x+ cos \- 1 which is small, of order 
Z 2 . Finally, the scattering phase 8x l = X+ ~ X- in (0) 
may be identified with Sx ti = X+ — X- j the deviation of 
5x l from 5x ti provides only small corrections of order Z 2 
to the splitting in the Andreev spectrum and thus can be 
disregarded. 
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FIG. 3. Discrete energy spectrum for a smooth parabolic 
potential barrier bounded by Dirac scatterers of strength Vq 
(Z = mVo/?i 2 fcp; parameters are chosen to emphasize the 
overall structure of the spectrum). In regions I and II the 
quasi-particle energies depend on the phase ip (ip = 0: solid 
lines, ip = 7r: dashed lines). Note that the degeneracies in 
region I are lifted for ip 7^ 0, ±7T. (a) The unperturbed spec- 
trum for Z — as discussed by Chthelkatchev et alx3 (b) At 
Z — 0.1 the junction exhibits weak resonances which lead to 
a splitting in the Andreev spectrum; the period A reflects the 
resonance spacing, (c) For Z = 1 the junction develops dis- 
tinct resonances to which the quasi-particle levels are pinned. 
The remaining degeneracies are a consequence of the almost 
identical resonance shapes. 
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At cp = the right hand side of (0) can be expressed 
as R = 1 - 2Z 2 (1 - cos5x**)(l - cos£x*0 while the 
left hand side takes the form L = cos(Sx ti — a)- The 
sum T>x li = X+ + X- 01 the scattering phases depends 
only weakly on e but increases linearly in fj, X Q with a 
slope 27r/A, as follows from linearizing x±( £ i l^xo)- On 
the other hand, the phase difference 6x li increases with 
energy e but is roughly independent of [i X Q. Both terms L 
and R then show an oscillating dependence on Sx** and 
on the energy e (see Fig. 4), while their relative phase 
a(e) slowly decreases from n at e = to at e = A. 
Thus, for small energies (e <C A) the maxima in L coin- 
cide with the minima in R and we obtain a large splitting, 

5e^-ZK\sin X %\, (8) 

7T 

to lowest order in Z (with x+ evaluated at e = 0). The 
corresponding splitting for ip = ±tt is proportional to 
|cosx+| with the same prefactor. With increasing en- 
ergy the phases of L and R match up and the splitting- 
becomes small, see Fig. 4. The monotonous decrease of 
the splitting is most prominent in long junctions with 
many trapped levels. 

As manifested in (|[) the splitting vanishes at specific 
degeneracy points: we define the (I, m)-degeneracy via 
the condition e F = (E l ICS + £™ s )/2, where I and m count 
the perfect resonances. This implies l-E'cs — = 
e ros = e rcs = ef cs . The ideal degeneracies in (|j) originate 
from the assumption that all perfect resonances have the 
same symmetric shape, tlE 1 ^ + SE) = t(E™ s ± 5E). For 
ip = (±7r) and for even (odd) values of n(e^ cs ) the right 
hand side of the quantization condition (multiplied by 
— 1) equals unity, r + r_ + t + t- = r\ + 1\ = 1, and thus 
the trapped levels appear in degenerate pairs. The de- 
generacy is lifted, if we take into account that the perfect 
resonances are not identical. 




A~" e 

FIG. 4. The quasi-particle energies are determined by the 
intersections of the right and the left hand side (R and L) 
of the quantization condition (^). Every maximum in L con- 
tributes a pair of trapped levels split by Se. 

B. Strong Resonances 

Let us next study the effects of strong resonances. 
Without loss of generality we again discuss the sym- 
metric barrier bounded by 5-scatterers (see also Ref. 
|23| ) . We concentrate on the open channels in region 
I (see Fig. 3(a)); the other regimes in the (e, fj, x o)-plane 
are less interesting and we will briefly discuss them at 



the end. For sharp resonances with V <C A the trans- 
mittivity close to a resonance e ros assumes a Lorentz 
profiled t±(e) = (T/2)/y/(e - e ICS ) 2 + (r/2) 2 and the 
scattering phase takes the usual form x±( £ ) = X*(-&res)± 
arctan[2(e — e rcs )/r]. 

Fig. 5 shows the terms entering the quantization con- 
dition ([?]). Particles incident from the left on the normal 
region are reflected back (r + sw 1) unless their energy co- 
incides with an electronic resonance energy (r + (e rcs ) = 
for perfect resonances). Similarly, holes are only trans- 
mitted if their energy corresponds to a hole-like reso- 
nance energy (r_(e rcs ) = 0). Consequently, the product 
r + r_ remains close to unity but sharply drops to zero 
at perfect resonance energies (see Fig. 5(a)). Similar to 
the case of weak resonances the product r + r_ becomes a 
smooth function of e when combined with cos (5 = (— 1)", 
see Fig. 5(b). Furthermore, the phase sensitive term 
t + t- cos ip vanishes practically over the entire interval 
[0, A] (see Fig. 5(c)) with the exception of specific points 
where particle- and hole- resonances become degenerate, 

^res — £res — ^res " 



(a) (b) (c) 




FIG. 5. The various terms entering the quantization con- 
dition (^) in the limit T <C A: We show the case of perfect 
isolated resonances (je — e\ 3> T), where the product t+t- is 
small at all energies (here, t+t- < 0.025). The bound state 
energies (marked with crosses in (e)) are determined by the 
intersections of cos(<5x' — a) and r+r_ cos/3. 

The left hand side of (0) is mainly determined by 
the phase difference 5% exhibiting sharp steps by 7r at 
the resonance energies e ros and e ros while remaining con- 
stant in between (in our discussion of the situation away 
from region I below, we will have to distinguish perfect 
from imperfect resonances, as the latter involve phase 
jumps by 27r). Hence, in the limit T — ► we obtain 
cos((5x* — cy) — (— l)"cosa while a small finite value of 
r will smooth the discontinuities (see Fig. 5(dV). All 
the terms entering the quantization condition (|7]) then 
are pronounced functions of e near the resonance ener- 
gies while staying roughly constant everywhere else. As 
a consequence, the intersections of the left and the right 
hand side of (m) come to lie close to the resonance en- 



5 



ergies of the NININ-j unction as shown in Fig. 5(e) and 
we conclude that the normal state resonances attract the 
quasi-particle bound states, e « e res to zeroth order in 
r. The predicate electronic- and hole-like states then can 
be naturally assigned to the trapped levels in the SINIS 
junction, too. 

In region I (see Fig. 3(a)), the bound state energy close 
to an isolated resonance is determined by the implicit 
equation 



• cot[a(e)/2] 



(9) 



This equation follows directly from the quantization 
condition (0) with t + exp(i%+) described by an ideal 
Lorentzian resonance while t_ = 0, r_ =1, and \_ — 
X^(e res ) + n(e rcs )Tr is constant (or vice versa + «-> — ). 
For e <C A we may approximate a(e) ~ tt, while at 
e < A (||) is easily analyzed graphically (note, however, 
that the above assumptions for t± and x± are only valid 
if |e — e res | <C A). The phase independence of the quasi- 
particle states in the presence of strong normal scattering 
is made explicit in the result (||). 

Close to the particle-hole degenerate resonances the 



phase sensitive term = (T/2) 2 /[(s — 



(T/2Y 



is of order unity and we obtain pairs of trapped levels 



r(±): 



£( co t(f)± 



£(cot(§)± 



sin(a/2) 
cos((/?/2) | 
sin(ct/2) 



1 1 n(ef cs ) even, 
! , odd 



(10) 



FIG. 6. Discrete spectrum for a triple barrier generating 
narrow resonances (from a numerical solution of the quan- 
tization condition (T7|)). The levels are pinned to the reso- 
nances; note that hole-like resonance energies £ res grow with 
(j, x o whereas electronic resonance energies e res decrease. In re- 
gion II (III) the hole-like (electronic) levels pair up such that 
the branches of imperfect resonances carry two bound states. 
The inset shows the lowest level near the degeneracy point 
Uxo f° r different values of f> (ip = j7r/4, j = . . . 4). 

Above we have found the trapped levels belonging to a 
branch £ r es(Aro) of perfect resonances at the degeneracy 
points fi^o and far away, \u X Q — M^ol ^- ^ n between we 
can interpolate (^) and ( |l0| ) by means of the usual hyper- 
bolic dispersion relation and obtain (in the limit £« A; 
a w 7r) 



,(±) 



(Mxo) 



e Q ± 

°rcs 



(^o-/4 ) 2 + [rcos(^/2)/2] 2 , (11) 



, rcs , is odd and d M:c0 e res 



where we have assumed that n(ef cs 

— <9 M[c0 £ r es ~ 1. The result ( [Tl"| ) will be useful in our dis- 
cussion of the transport properties below. 

In region II (see Fig. 3(a)) the perfect hole-like reso- 
nances pair up and collapse as /z^o is lowered, thus gener- 
ating imperfect resonances, see Fig. 6. The analogous col- 
lapse of the electronic resonances is shifted to region III. 
Note that an imperfect resonance carries two nearly de- 
generate levels guaranteeing that the number of bound 
states remains conserved upon changing [j, x q. 



C. Other Sources of Non-ideal Transmission 



(we assume both scattering amplitudes t± exp(«Xj_) to be 
described by the same Lorentzian centered around ef cs ; 
Eq. [l(] aspees with the result of Ref . || obtained for a short 
junctionEZl, see also Ref. ||). The levels (|l(|) are mani- 
festly phase sensitive and become degenerate at <p = for 
even n and at tp = it if n is odd. For low energies e <C A 
we estimate a(s) ~ tt and the maximum level splitting 
is Se » T. We conclude that phase sensitivity survives 
only in a narrow interval of order A[i x q ~ T around the 
degeneracy points /j, x0 , being negligible everywhere else 
(see Fig. 3(c)). 




The above analysis has been based on <5-scatterers 
modelling the effects of an insulating layer in an SIN 
interface. Assuming other sources of non-ideal trans- 
mission at ±L/2 we have to modify some of the above 
results: In order to describe junctions with potential 
steps V s at the edges we introduce the wave vector 
k nw = yj2m(E — Vgj/h in the normal wire and distin- 
guish it from the analogous quantity fc sc = V2mE/h in 
the superconductor. With n = (k sc / k nw + k nw / k sc ) / 2 > 1 
the transmission probability takes the form T = 2/(1 + k) 
(see the Appendix) . In the presence of weak resonances 
(T < 1), Eq. (||) for the gaps in the Andreev spectrum 
reads Se ps AAy/T^T\ cosx**!/""; while the results for 
strong resonances remain the same (see (0)). A step in 
the effective mass at the interface produces a similar re- 
sult: with m(\x\ > L/2) = m sc and m(\x\ < L/2) = m nw 
the ratio of wave vectors entering the parameter k takes 
the form k sc /k nw — y/ m sc /m nw . Finally, combining all 
three effects, a potential step, an effective mass disconti- 
nuity, and an insulating layer (in the form of a 5-function 
scatterer) at the interface, the transmission probability 
can be written in the form 



(12) 



m nw ) 2 + W 2 ml c m1 w /Ti 



G 



2min ^ Tmax, exhibits large imperfect 



obviously, the various scattering mechanisms are non- 
additive and Matthiessen's rule is not applicable. 

Slightly asymmetric junctions have essentially the 
same properties as the idealized symmetric ones. An 
NINFN junction with insulating layers of different 
transparei 
resonance: 



of height t 2 {E rcs ) « T min /r max instead of 
perfect resonances with unit transmission. Furthermore, 
X r is now a continuous function of energy; the discontin- 
uous jumps by tt are smeared (with \ r increasing mono- 
tonuously for T_ L / 2 < T L / 2 and bounded for the opposite 
case with T L / 2 < T_ L / 2 ; note that \ r + X r = n + 2%*, 
X r the corresponding phase for a particle incident from 
the left) and cos/3 changes rapidly but smoothly be- 
tween ±1 at the resonance energies. Since all terms 
entering the quantization condition (Q) are only weakly 
affected by small asymmetries, the above discussion of 
the electronic properties remains valid for (weakly) non- 
symmetric junctions. 



IV. TRANSPORT 

We proceed with the investigation of the transport 
properties of the symmetric SINIS junctions in region 
I (see Fig. 6) and concentrate on the situation charac- 
terized by strong resonances with T/A <C 1. Given the 
dependence of the quasi-particle energy e on the phase if, 
the contribution of the level to the supercurrent follows 
from a simple derivative, / = (2e/h)d v e (a factor 2 has 
been included to account for spin-degeneracy). 



A. Generic case 



contributes with a term of order T 2 . The only situa- 
tion producing a large current (of order T) then is re- 
alized at the special degeneracy points produced by an 
(electronic) resonance crossing the Fermi level, where 
£ res (/4o) « e res (Mxo) ~ °' see Fi g- 6 - Witri its energy 
£q + ^ described by (|ll|), this level carries a nonvanishing 
supercurrent of magnitude 



la = 



2eT 2 



sin tp 



!6 v'(^o-^o) 2 + (r/2) 2 cos2(^/2) 



(13) 



With one channel open, we then find that a large critical 
supercurrent 



eT 
2h 



Ox0-/4o)l 2 iMzO 



r/2 



Mxol 



r/2 



(14) 



is realized near the special values /j, x0 for the chemi- 
cal potential where sf es {^xo) = and for ip = n — 0. 
This large supercurrent flow quickly vanishes as /i x o is 
tuned away from these degeneracy points by an energy 
larger than the resonance width T, see Fig. 7. Below, 
we refer to such a tunable SINIS junction as a Resonant 
Josephson-Transistor (RJT) (alternative schemes leading 
to a transistor-effect make use of resonant electromag- 
netic pumpingE3 or injection of quasi-particlcs into spe- 
cific levels via multiprobe devices, see Refs. |28|,|30|,[l5] fL7| ). 
Note that the chemical potential (j, x q is not directly ac- 
cessible but only through the gate voltage V g , thus in- 
troducing the slope d[i X Q/dV g as an additional character- 
istic parameter of the device. Also, we point out that 
a sharp transistor effect requires the temperature to be 
low, k B T < r. 



We first consider the generic case where a level is 
pinned to an isolated (electronic) resonance. In this situ- 
ation we can assume t+ ~ 1, t_ oc T/A and cos(^x*~ a )~ 
r + r- cos/3 oc e/r + const, (see also Fig. 5). Within this 
approximation oc T/A is small and linearly related 
to Se — max v [e((/))] — min^fe^)] through the large slope 
1/r. We then estimate Se oc d v e oc T 2 and find that each 
level contributes a small supercurrent, of order T 2 . 



B. Degenerate resonances 

The phase sensitivity of the trapped levels is dramat- 
ically increased close to degenerate resonances, see Eq. 
(|lT|). However, still assuming narrow resonances with 
T/A <C 1, the supercurrent in general remains small: 
within our approximation the contributions arising from 
a pair of nearly degenerate levels cancel each other due to 
the symmetry d v e^ = —d v e^' (a more accurate anal- 



ysis provides a residual contribution of order T , seeEf). 
Furthermore, the continuous part of the spectrum again 
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FIG. 7. Critical current versus chemical potential 
Uxo = 8h/2r\ £F for a sharp resonance with T = 0.1 and 
t| £f = Wh/A. The approximation dli| ) (dashed line) agrees 
with the exact result obtained via the Green's function anal- 
ysis at resonance and midway between resonances, 8 = 0, ±7r, 
and provides a good estimate in between. 
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The result (g|) appears to be valid in the immediate 
vicinity of a resonance crossing the Fermi level. How- 
ever, comparing this result with- the one derived from 
a full Green's function analysisEJ of the long [L S> £) 
SINIS junction with sharp resonances (r/A <C 1) one 
finds that the result (|l4j ) is exact at the degeneracy 
points HxQ ~ Uxo an( i, as it turns out, also midway in 
between two such degeneracies; defining the parameter 
6 = 2S\ 6F /H w 2(/x x0 - fJ-io)T\ ep /h + 2mr the result (p|) 
then is exact at all multiples of it, S sa n7r (here, 5[7 F 
denotes the action for an electron traversing the normal 
region [— L/2 + 0, L/2 — 0] at the Fermi level and r = d e S 
is the travel time). Furthermore, in between the points 
5 = mr the deviations of (Q) from the exact Green's 
function result are small, see Fig. 7. 



C. Josephson- versus Single Electron Transistor 

Next we discuss the relation between the Josephson- 
Transistor and the Superconducting Single Electron 
Transistor (SSET)l23. The latter consists of two Joseph- 
son junctions separated by a superconducting grain, usu- 
ally referred to as the island, which is capacitively (with 
capacitance C g ) coupled to an external gate electrode. In 
a more figurative terminology we might call this device a 
SISIS junction. While the charging energy Ec = e 2 /2Cs 
tends to fix the number of Cooper pairs on the island, 
the (conjugate) phase variable tends to be fixed by the 
Josephson coupling energy Ej of the two junctions. For 
the SSET, charging effects are dominant, E c ^> Ej, and 
thus the total capacitance Cs of the island must be small 
(but stilEa C g <C Cs). Moreover we assume Ec < A 
and hence the ground state of the island contains an 
even number of electrons at any gate voltage. Given this 
layout for the device, we expect the supercurrent to be 
suppressed by the 'Coulomb blockade' effectEJ. However, 
for specific values of the gate voltage the energies of 
two even charge states differing by 2e become degener- 
ate. At these 'Coulomb resonances' the suppression of 
the Josephson current through the island is lifted and 
we find a large superflow. We can control the supercur- 
rent by changing the gate voltage and hence the SSET 
constitutes a transistor device. 

The above discussion shows that the SSET and the 
RJT are like devices, switching on as levels become de- 
generate; while for the SSET these levels belong to fixed 
charge states of the island, for the RJT these levels de- 
rive from transmission resonances. In both cases the de- 
generacy is lifted by the Josephson coupling E cos(<p/2), 
where the coupling constant E is related to the trans- 
mission of the NIS boundaries, see below. Furthermore, 
the SSET and the RJT exhibit largely the same current- 
voltage characteristic close to these level degeneracies, 



2e E 2 



sm tp 



h 4e7 



yJ(y g -V*)* + (Ecos(<p/2)/ejf 



(15) 



where the parameters E and 7 have to be specified for 
each device. In fast, (15) yields the current through 
a symmetric SSET23 if we substitute the energy pa- 
rameter E by the coupling energy Ej and the dimen- 
sionless constant 7 by the small xatio Cg/Cg. Using 
the Ambegaokar-Baratoff relationtij and the Landauer 
formulae^ we can reexpress the coupling energy in terms 
of more microscopic quantities: Ej = ATA/4, where 
N (x k 2 A denotes the number of open channels in the 
tunnel junction, A denotes its area, and T its transmis- 
sion coefficient. On the other hand, the RJT involves 
the parameters E = T/2 and 7 = dfi x o/deV g . Again, 
the resonance width is determined by geometrical quan- 
tities: r w TA/ir, where T is thej-transmittivity of the 
insulator layer in the triple barriertll. For flat inner bar- 
riers the quasi-classical method provides the simplifica- 
tion A ~ 2(fi x0 e Ij ) 1/2 , with e L = h 2 Tr 2 /2mL 2 . Note that, 
while charging effects dominate the physics of the SSET 
these are much less relevant for the RJT as the latter in- 
volves an 'open wire' rather than a 'closed island'; the 
adiabatic joints to the superconducting banks provide 
reservoirs which effectively screen the charge transport 
through the wire. 

From Eq. ( |l5| ) we can determine the critical supercur- 
rent I c as a function of the gate voltage near V«. Appar- 
ently, I c (V g ) defines a current peak of width 5V g = E/ej 
that attains its maximum 7 max = eE/% at the degener- 
acy point. We estimate its slope by the ratio I ma x/5V & = 
e 2r y/h carrying the dimension of a conductance. For the 
SSET 7 -c 1 is determined through the small capaci- 
tance ratio Cg/Cg. The analogous quantity for the RJT 
is of order unity, implying that the RJT shows a more 
prominent slope: with dfj, x o/deV s = dd/^xo / 'd e v„d, and 
ddl^xo = 2(e F — Hxo)/d we have to determine the sucep- 
tibility of the channel width d with respect to the gate 
voltage V g . Making use of simple electrostatic consider- 
ations one easily finds that d^ d ~ a B /e F , where a p 
notes the Bohr radius in the semiconductor materialHI 
Hence, 7 ~ a B /d and with a B of order 10 nm typically 
we arrive at a value of order unity for the parameter 7. 

Comparable energy scales A ~ \/JkM£Z can be 
reached in realistic short SINIS junctions (L < 
Tr(lJ>xo/2m) 1 / 2 h/ A). Furthermore, we can optimize the 
performance of the RJT by tuning the transparency of 
the insulating layers (the Dirac <5-scatterers) until the 
resonance width T approaches the resonance spacing A; 
such a choice of parameters produces still isolated current 
peaks of maximum height. On the other hand, the maxi- 
mum superflow through the SSET scales with the number 
N of open channels, I ma x = eNTA/ATi, while the RJT as 
defined above is generically a single channel device with a 
critical supercurrent 7 max = eTy/fi x o£ L /hir; going over to 
a many channel RJT device, resonances from individual 
channels superpose at different values for the_chemical 
potential and hence do not add up in generalo. Also, a 
smoothing of the resonance structure has-,been observed 
in the numerical results by Wendin et alB 
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Summarizing, superconducting transistors can be de- 
signed in terms of 'charge' (the SSET) or 'phase' devices 
(the JT). In the charge device, the island is separated 
from the superconducting leads through insulating bar- 
riers. On the other hand, switching channels in a perfect 
SNS junction one obtains a phase device with a (large) 
critical current I c — e/(ro + fij A) determined by the ef- 
fective time,-Sf) + h/A the charge needs to traverse the 
normal wiretJ. Such devices with perfect interfaces are 
difficult to fabricate — in reality we always have to ac- 
count for a non-ideal transmission T < 1 through the 
SIN interfaces and we end up with a Fabry-Perot type 
resonator device. Sequential tunneling through the in- 
terface barriers then reduces the supercurrent by a factor 
T 2 in general. However, tuning the chemical potential to 
a scattering resonance the critical current of the junction 
remains large, I c ~ eT/2h of order T — again, the crit- 
ical current is given by the time h/T the charge spends 
in the junction. 

A second result we wish to emphasize here concerns 
the fact that the critical current is carried by the lowest 
quasi-particle level alone. This has been demonstrated 
for the SNS device in Ref. and above for the SI- 
NIS junction for the case of strong resonances (devia- 
tions from the single level result are numerically (but not 
parametrically) larger away from weak resonances). The 
importance to know the dispersion of this level then pro- 
vides a good a posteriori reason for studying the quasi- 
particle spectrum in such junctions. Corresponding spec- 
troscopic experiments can be realized using multiprobe 
devices as proposed in the work of van Wees et allEi 
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APPENDIX: RESONANCE STRUCTURES IN 
ONE DIMENSION 



We consider a one-dimensional potential landscape 
V(x) confined to the interval [-L/2, L/2]. As a first 
example we concentrate on a symmetric barrier con- 
sisting of a broad (mfl 2 = -d%V, frQ < V(0)) and 
smooth potential barrier, bounded by steps of height 
V a = V(L/2-0)-V(L/2+Q) (see the inset of Fig. 8). Be- 
low we make use of the wave vectors k s , 



\ / 2mE/h and 

y/2m(E — V s )/h describing particles of energy E 
in the wire close to the boundaries. The global trans- 
mission amplitude iexp(ix*) iSr«ost easily obtained by 
determining the transfer-matrixE-J of the barrier between 
—L/2 — and L/2 + and we find the result 

texp{i X t ) = rT r . -. , t ... , (Al) 

cos(x *) — iKsm(x ') ± W"i<7 



where tjexp^x**) and riexp(ix Vi ) denote the transmis- 
sion and reflection amplitudes of the inner smooth bar- 
rier. In (Al) we have introduced the definitions k = 
{k sc /k DW + k nw /k sc )/2 and o~ = (fec/tw k nv/ /k sc )/2 
containing the information about the steps. In our ex- 
ample, k, a, and U show a weak dependence on energy 
and remain almost constan t w hile exp(ix* i ) oscillates 
rapidly. The derivation of (Al) makes use of the rela- 



tion x 



X' 



tt/2 



n G Z, which follows from 



the unitarity of the scattering matrix; the term ±ir^a 
changes sign at each perfect resonance of ti. Since the 
inner barrier was assumed to be smooth, we expect ti 
to exhibit no resonances at all and the term r^cr always 
carries a positive sign. 



Next, we show that (Al) qualitatively reproduces the 
transmission amplitude shown in Fig. 2. We definine the 
function z — l/iexp(zx*), whose four arguments x ■> k j 
<t, and ti depend on the energy. Since the inner scat- 
tering phase x li is the most energy sensitive argument 
we minimize \z\ with respect to x ti (keeping k, er, and 
ti fixed at a given energy E). Thus we obtain the sta- 
tionary phases Xmin(^) belonging to the minima of \z\. 
Moreover, we can introduce the complex valued func- 
tion z miu [E] = z[x^ in (E),K(E),a(E),U(E)], which is 
roughly constant when compared to z[E]. We can es- 
timate the resonance energies by solving the equation 
z[-E res ] = •Zmin[-E'rcs] for E rcs . Apparently z[E] describes 
an ellipse in the complex plane with half axes 1/ti and 
n/ti. The center of this ellipse is shifted away from the 
origin by iria/U. 



Im(z) 




' Re(z) 







IhQ. r- 




J s 





-U2 


L/2 x 



h < fcri. < 1 



FIG. 8. Complex transmission amplitude through a 
smooth barrier bounded by steps. With decreasing energy 
E the inverse z = l/£exp(ix 4 ) moves counterclockwise on an 
ellipse, the center of which is positioned at ina/ti and shifts 
up the imaginary axis. As ti drops below t cr i t the perfect 
resonances at P and Q collapse at the bottom of the ellipse 
(P — > S <— Q) and transform into an imperfect (double) res- 
onance. 

Fig. 8 illustrates the behavior of z m ; n when grows 
from to 1 in the intermediate region E E [V(0) — 
h£l/2, V(0) + hfl/2] (with a smooth inner barrier, is 
strictly monotonous, r.i ps for E > V(0) + hfl/2 and 







n w 1 for E < V(0) - The distance from the 

origin to the bottom of the ellipse (point S) is always 
extremal. At high energies where 1 ~ i; > l/n = t cr it, 
we observe perfect resonances: with |z m i n | = 1 these res- 
onances are realized at the symmetric points P and Q 
where the ellipse touches the unit circle, see Fig. 8. As U 
decreases the perfect resonances approach each other (P 
and Q move towards S) and merge with S as U^t CT it. At 
energies E < V(0) (=$■ U < t CI a) the reflection coefficient 
Ti drops to zero, and S becomes the closest point to the 
origin with |z m i n | = (/t— 0Tj)/t, > 1; the resonances then 
have paired up and have become imperfect, with a height 
decaying rapidly with decreasing energy E. 

Physically, this resonance structure originates from the 
interplay of two competing transparencies. The trans- 
mission probability through the potential steps is almost 
constant and much smaller than the transparency of the 
inner barrier for E 3> V(0). In this regime the parti- 
cles propagate freely and are reflected by the steps alone; 
thus we observe perfect resonances, and the shape of the 
smooth inner barrier only influences the resonance spac- 
ing. On the other hand, for E <C V(0) the inner barrier 
behaves like a hard wall while the steps are compara- 
tively transparent. The global transmission through the 
barrier bounded by the steps is suppressed by the tun- 
neling through the classically inaccessible region. The 
resonances of the subsystems consisting of a step and a 
potential hill account for the global imperfect resonances. 
Note that the barrier bounded by steps generates pro- 
nounced resonances only if V s > 0.7 V(0). 
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FIG. 9. The solid curve in the (r ,n) plane characterizes 
the energy dependence of the inner and the global reflection 
amplitudes in a symmetric triple barrier as shown in the in- 
set. The domains of perfect and imperfect resonances are 
separated by the curve r cr i t (r ). In an narrow interval ~ Ml 
around V(0) the reflection coefficent r 2 drops from 1 to 0; in 
this energy range r can be considered as a constant. 

In a second step we investigate the transmission am- 
plitude through three arbitrary symmetric potential bar- 
riers, i.e., a triple barrier, see Fig. 9. The inner barrier 
is again smooth, whereas the two outer barriers (char- 



acterized by the transmission and reflection amplitudes 
t exp(ix to ) and r exp(ix r °)) are assumed to be equal. 
Thus the global transmission amplitude takes the form 



te lx = 



1 + 2r. l r e i (x t *+x t °) - r 2 el 2( x *.+ x 'o) • 



(A2) 



In order to determine the resonance structure we eval- 
uate the extrema of the modulus of the denominator in 
flA2] ) with respect to \ li + ( we nx h an d t a , which 
are less sensitive to changes in energy) . A direct calcula- 
tion yields the two extremal conditions cos(x* i + X*" ) = 
-n(l + r 2 )/2r (I) and sin(x 4, + x*°) = ( n )- Elemen- 
tary manipulations show that the transmission amplitude 
exhibits perfect resonances in the first case. But condi- 
tion I can only be satisfied as long as rt < r cr it(r ) = 
2r Q /(l +r~), i.e., for a reasonably transparent inner bar- 
rier (see Fig. 9) . The extremal condition II requires a sec- 



ond distinction. Condition 11(a) reads \ 



X 



27m, 



n G Z. The critical points of this type always belong 
to local minima of t. Finally, we discuss condition 11(b), 



it + 27m, n £ Z, which refers to local minima 



when rj < r cr it and characterizes imperfect resonances 
in the regime Ti > r cr it- Since the internal barrier was 
assumed to be smooth, the regime of perfect resonances 
lies at high energies E > V(0) + M7/2, where « 
and r ^> r,. At low energies E < V(0) — HQ/2, the 
resonances are always imperfect. Within the crossover 
region the perfect resonances attract each other pairwise 
and collapse to become imperfect as the energy decreases 
(see Fig. 2). 

Other sources of resonances are discontinuities of the 
force d x V and effective mass steps. Technically they can 
be treated like Dirac <5-scatterers and potential steps, re- 
spectively, as can be easily checked making use of the 
transfer-matrix formalism. 
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